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ANALISE MATEMÁTICA IV 

FICHA SUPLEMENTAR 5 

EQUAÇÕES DIFERENCIAIS PARCIAIS 
E TRANSFORMADA DE LAPLACE 



Séries de Fourier 



(1) 



Desenvolva em 


série de Fourier as seguintes funções: 


(a) /:[-!,!] 


— > IR definida por 






/<*) = {;; 


se - 1 < x < 
se < x < 1 


(b) /:[-7T,7T] 


— > IR definida por f(x) 


= cos 2 (x) — sin 3 (x). 



Resolução: 

(a) A série de Fourier da função f é da forma 



«o 



+ y, (a n cos(mTx) + b n sm(mrx)) 



n=\ 



onde 



f(x) cos(n7ra;) dx 



para n > 



b n = / f(x) sin(n7rx) dx para n > 1 

Como f é uma função ímpar, todos os a n s são zero e 



b n = 2 I f(x) sin(n7rx) dx 



= 2 sin(n7ra;) dx 



— — (cosnir 

717T v 


— cos 0) 


2-2(-l) 7t 




ri7r 




í- 


se n é ímpar , 


o 


se n é par . 



Logo, a série de Fourier de f é 



E4 r— v 4 
— smímrx) = y — sin((2n+ lWx) . 



n = l 

n ímpar 



UTT 



- (2n + IH 



n=0 
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Comentário: Como a função / é seccionalmente diferenciável em [—1,1], a sua 
série de Fourier converge em [—1, 1] para 



f(x) se x g] - 1,0[ ou x g]0, 1[ 







se x = ou x = ±1 



-1 se x g] - 1,0[ 

se x = ou x = ±1 

1 se a; G]0, 1[ 



Nos pontos da forma x = k G Z, a série converge para porque este valor é a média 
dos limites laterais nestes pontos da extensão de / a R como função periódica de 
período 2: 



lim f(x) = lim f(x) = 1 e lim /(a?) = lim f(x) = — 1 . 

x — > — 1~ x — >1 — a: — >1+ x — > — 1 + 



(b) A série de Fourier da função f é da forma 







«o 



+ /^ (a ra cos(nx) + ò ra sin(na;)) 



ra=l 



L/ma vez que o desenvolvimento de Fourier é único, qualquer desenvolvimento que 
obtenhamos para f em termos das funções base 1, cos(nx) e sin(nx) será o desen- 
volvimento de Fourier de f . Ora tem-se 

2 1 + cos(2a;) 11 

cos x = = — I — cos(2a;) 

2 2 2 v ' 

e, pela fórmula de DeMoivre, 

Im ((cosx + isinx) ) = Im(cos(3x) + ísin(3a;)) 
<í=> 3 cos a; sina; — sin x = sin(3ai) 

<í=^> 3(1 — sin x) sin x — sin x = sin(3x) 

^ 3 1 . , 

•^=^> sm a; = -sina; — -sm(3x) 
4 4 

portanto o desenvolvimento de Fourier de f é 

11..3 1 , , 

- + - cos(2a;) — - sin x + - sin(3a;) 

isto é, ao = 1, ã2 = \i a n = para n ^ 0,2 e b\ = — |, 63 = ^, ò n = para 
n^l,3. 

D 



(2) 



Desenvolva 


as segu 


intes funções em série 


de Fourier: 


(a) m - 


í e x 

~ \ e~ x 


se < a; < 1 

se - 1 < x < ' 




(b) / : [- 


7T,7r] -» 


IR definida por f(x) = 


sin 4 x ; 


(c) m - 


í X 

10 


se - 2 < z < 
se < x < 2 
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Desenvolva a função / 


[0,7T]- 


-> IR definida por 

f(x) = e-* 


em série de senos. 







(3) 

Resolução: A série de senos da função f é da forma 

oo 

^ K sin(? 
onde, para n > 1 se tem 



[nx) 

re=l 



b n = 


2 r 


/(a;) sin(nx) da; 




2 /■"■ 


e~ x sin(nx) dx 




TT 7o 






— Im 


í r ei _ 1+in)x dx 




TT 


\Jo 




— Im 

7T 


/ e (-l+ire)a; 


'1 




l — 1 + ín 


„J 




— Im 

7T 


fe _,r (-l) n -l\ 




V -1 + m / 




2n(l 


-e"^(-l) n ) 



7r(l + n 2 ) 
/.ogo, a sér/e de senos de f é 

^ 2n(l-e-*(-ir) . 

> t^ ^ sm(nx). 

^ 7r(l + n 2 ) 



D 



(4) 



Seja l um número real positivo. Desenvolva a função / : [0, l] — > K. definida por 

x se < x < ; 

2 

em série de cosenos. 



1 i - x se £ < x < l 



Resolução: A série de cosenos da função f é da forma 

oo 

_ + XX 

onde, para n > se tem 



2 

ra=l 



oo 



( '" ' x ' cos'""'" 



2 /" z „, N n7rx 



fln = y / /(a:) cos — - — da: 



/4ss;'m, 



1 -z 



2 / /"2 /■' \ l 

ao = - / a; dx + / (/ — x) dx I = - 
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e, para n > 1, 



a„ 



= jl í 2 xcos(^f-) dx+ í(l-x)cos(^f-) dx 
2 Ix " ' ' 



-siní^ 
n7r y 



717T 



./o "7T 

ri 



sin(=p) da; 



/ -— sin(^) <fe 
/i n7r 



I 2 / / 2 \ 

- S in(^)- cos(=p) 

2n7r 2 \7r7H J / 



2/ 

= ^2( c °s(f)-l-co S (n 7 r) + co S (f)) 

2/ 
= -^(2cos(^)-l-(-]^ 

Tendo em conta que 



cos (f ) 



tem-se 



se n é ímpar 

1 se n = 4fc para algum ii 

-1 sen = 4fc + 2 

— |^ sen = 4/c + 2 




inteiro k 



— fr se n = 4/c + z 

CL = \ "" 

\ caso contrário 

Portanto a série de cosenos de f é 

]_ , V~^ 8/ f (4fc+2)7ra \ 

4 + Z. (4fc + 2)V C °H i ; 



(5) 



(6) 



D 



Desenvolva 


as segu 


intes fur 


ções 


em 


série de 


cosenos: 


(a) m = 


■{i- 


se 
2x se 


< 


x < 
x < 


í 

2 para 


x G [0,1]; 


(b) f(x) = 


= cos 2 X 


— sin 2 x 


para 


x G 


[o,f]. 





Desenvolva 


as seguintes funções em série 


de 


senos: 


(a) f(x) = 


= 1 pa 


ra x G [0, 1]; 








(b) m = 


= (T" 


x)x para x 


G[0,7r]; 






(c) m = 


= 55- 


40x para x 


e[0,l]; 






(d) f(x) = 


= x(x 2 


— 1) para x 


e[0,l]. 







Comentário: /\s respostas às primeiras três alíneas deste exercício estão contidas nas 
resoluções dos exercícios 1, 13 e 9 respectivamente. <) 
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Equações do calor, de Laplace e das ondas 



Seja c um parâmetro real positivo. Recorrendo ao método de separação de variáveis, 
determine a solução do seguinte problema para a equação das ondas 



(7) 



d 2 u _ c 2d 2 u 
dt 2 J dx 2 



u(t,0) = u(t,l) = 
u(0,x) = , Íf(0,x) = 1 

para í>0e para x e [0, 1] (satisfazendo a equação diferencial para < x < 1). 



Resolução: Pelo método de separação de variáveis, procura-se soluções da EDP da forma 

u(t,x) = T{t)X{x) , 
para as quais 

di 2 



d \T(t)X(x) = T^(t)X(x) 



£,T(t)X(x) = T(t)X^(x). 
Substituindo na equação, e assumindo que T(t)X(x) ^ 0, fica 

T^(t)X(x) = c^T(t)X^(x) <=► ^| 1 = C 2 ^f 

k k 

onde cada um dos membros tem que ser uma constante real k (porque o membro esquerdo 
não depende de x, o membro direito não depende dete eles têm que ser iguais). 
Para que a condição na fronteira, 

T(t)X{0) = T(t)X{l) = , 

seja satisfeita por uma solução não identicamente nula, tem que ser 

X(0) = X(l) = o 

(já que, se X(0) ^ ou X(l) ^ 0, teria que ser T(t) = 0, Víj. 



A solução geral de 

~2 



X {2 \x) = \x(x) 
c z 



é, consoante os casos, 



X(x) = c\e~ x + c 2 e~~ x se k > 

X(x) = c\ + C2X se k = 

,X(x) = cicos(^^x) + c 2 sin(^=^a;) se k < . 

Quando k > 0, a única solução X{x) que satisfaz a condição na fronteira é a solução 
identicamente nula, ou seja, nesses casos, 

X(0) = X(l) = => d = c 2 = . 

Quando k < 0, a condição X(0) = impõe c\ = e depois a condição X(l) = impõe 
a uma solução não identicamente nula que 

sin- ) = . 

c 

Esta condição pode ser satisfeita por uma solução não identicamente nula desde que 

\f—h 

= nn , n=l,2,3,... 



ou seja, 
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k = —c n 7T , n= 1,2,3,... 



Para k < O, a solução geral de 



r^(í) = JfeT(í) 

é 

T(t) = c\ cos(V-kt) + C2 sin(v / — kt) . 

A condição inicial 

T(0)X(x) = O 

impõe que seja T(O) = O para que a solução não seja identicamente nula. Quando T{t) é 
uma combinação linear de sin e cos como acima, tem que então ser c\ = 0. 

Obtêm-se assim todas as seguintes soluções para o problema de valor na fronteira con- 
jugado com a primeira condição inicial: 

u n (t, x) = sm(cmrt) sin(n7nc) , n = 1, 2, 3, . . . 

T(t) X{x) 

definidas para t > e para x e [0,1]. Por linearidade, qualquer combinação linear finita 
desta soluções é ainda solução e, mais geralmente, qualquer série 



/ "n^nK^i '• 



n=\ 



é uma solução formal. 

Para que a solução mais geral da forma acima satisfaça a segunda condição inicial 



tem que ser 



Como 



a condição fica 



du 






(0,1) = 1 , 



du. 



dt 



= í . 



í=0 



du n 
~dt 



= cn-K cos(cmrí) sm(nirx) 



/ dnCn-K sm{n-Kx) = 1 . 



71=1 

Para encontrar as constantes d n adequadas, desenvolve-se a função 1 em série de se- 
nos. Isso é equivalente a estender esta função ao intervalo [— 1, 1] como função ímpar, e 
desenvolver a extensão em série de Fourier, trabalho esse que foi feito no exercício 1 onde 
se encontrou que a série de senos para a função 1 é 

oo . oo . 

E4 ^— v 4 
— sin(n7rx) = > -. : — sin((2n+ lWa;) . 

n = l n=0 v ' 

ri ímpar 

Logo, os d n 's devem ser quando n é par e devem satisfazer 

4 



d„cnir = 



titt 
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quando n é ímpar. 

Conclui-se que a série que dá a solução formal do problema posto é 



4 



— k-^7 sm(crnrt) sm(mrx) 



71 = 1 

n ímpar 



= / ~F7. , -» o o sin (c(2n + lWt) sin((2ra + l)irx) . 

^— ' c(2n + 1) Z 7T Z 

n=0 



D 



(8) 



Use o método de separação de variáveis para determinar a solução do seguinte 
problema de valor na fronteira para a equação de Laplace: 

= 



d 2 u _|_ d 2 u 
dx 2 dy 2 

u(x, 0) = u(x, 7r) = u(0, y) = 
u(ir,y) = sin(2y) 



para 0<x<7re0<y<7r. 



Resolução: Pelo método de separação de variáveis, procura-se soluções do problema da 
forma 

u(x,y) = X(x)Y(y). 
Substituindo na equação obtém-se, para X(x)Y(y) ^ 0, 

*»(,)n,) + *(,)y«(,) = o ~ ^ = -^f 

Conclui-se que existe k eK ia/ que 

x( 2 )(x) _ y( 2 )(y) 

Para qive a condição na fronteira 

X{x)Y(0) = X(x)Y(tt) = 

se/a satisfeita por uma solução não identicamente nula, tem que ser 

Y(0) = Y(n) = 

A equação 

Y^(y) + kY(y) = 

só tem soluções não identicamente nulas que verifiquem a condição fronteira Y(0) = 
Y(ir) = se k > 0. A/esse caso, a solução geral é 

Y(y) = ci cos(Vky) + c 2 sin(v ky) 

A condição Y(0) = implica c\ = e a condição Y(ir) = impõe 

sin(vK7r) = <í=> k = n 2 para algum n inteiro 

Para que a solução não seja identicamente nula, tem que ser n ^ e claramente basta 
considerar n > 0. A solução geral da equação 

X (2) (x)-n 2 X(x) = 

é 

X(x) = Cl e nx + c 2 e- nx 
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Para uma solução não identicamente nula, a condição na fronteira X(0)Y(y) = O implica 

X(0) = ^^ Cl + C 2 = O ^^ C2 = -ci 

pelo que 

X(x) = 2ci I J = 2ci sinh(nx) 

Obtêm-se assim as seguintes soluções que verificam todas as condições na fronteira excepto 
possivelmente a que foi imposta para x = n: 

u n (x, y) = sinh(nx) sin(ny), n = 1, 2, 3, . . . 



definidas para 0<x<ireO<y<Tr. 
Por linearidade obtem-se a solução formal 



y~]d n u n (x,y) 

71=1 



Para que esta série satisfaça a condição fronteira que resta devemos ter 



y^d n u n (ir,y) = sin2j/ 

ra=l 

oo 

y d n sinh(n7r) sin(ny) = sin 2y 



ra=l 

Conclui-se que d n =0 para n ^ 2 e que 

1 
£?2 = — 



sinh(27r) 
/sío é, a solução do problema do enunciado é 



U ( X >V) = ■ -um \ sinh(2a:)sin(2j/). 
smn(27rj 



D 
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Considere a equação do calor 



du r,d 2 u 

- a 1 



(9) 



dt dx 2 

onde a é um parâmetro real positivo, í>0eie [0, L]. Uma solução estacionária 
da equação do calor é uma solução que não depende do tempo, t. 

(a) Mostre que todas as soluções estacionárias da equação do calor são lineares 
em x, i.e., são funções da forma 

u(x) = ax + b . 

(b) Determine uma solução estacionária para a equação do calor com a seguinte 
condição na fronteira: 

u(t,0) = T 1 e u(t,L) = T 2 . 

(c) Ache uma solução da equação do calor (verificando a equação diferencial para 
< x < L) com as seguintes condições de fronteira e inicial 

ií(0, x) = 75 

u(t,0) = 20 , u(t, 1) = 60 . 

Sugestão: Considere soluções da forma 

u(t, x) = v(x) + w(t, x) , 

onde v(x) é uma solução estacionária do problema com condição na 
fronteira: 

u(t,0) = 20 e tt(í,l) = 60 . 



Resolução: 

(a) Seja u(x) uma solução estacionária da equação do calor. Substituindo u{x) na 
equação, obtém-se 

= a V 2 ) 
cuja solução geral é 

u(x) = Cl + ClX , 

que é uma função linear. 

(b) Procura-se os coeficientes c\ e c% tais que a correspondente solução estacionária, 

u(x) = ci + C2X, satisfaça 

fu(0) = ci = Ti 

\u(L) = ci + c 2 L = T 2 . 

A solução deste sistema linear é 

íci = Ti 
\c 2 =^. 

Encontra-se então a seguinte solução estacionária para a equação do calor com a 
condição na fronteira dada: 

í \ rp . T 2 - Tl 

U{X) = l\-\ X . 



(c) De acordo com a sugestão, procura-se uma solução da forma 

u(t, x) = v(x) + w(t, x) , 
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onde v(x) é uma solução estacionária do problema da alínea (b) tomando T\ = 20, 
T 2 = 60 e L = 1. Então 

v(x) = 20 + 40x , 

e w(t, x) = u(t, x) — v(x) terá de ser uma solução do problema homogéneo 

Í dw 2 d 2 w 
dt ~ a dx 2 
w(t,0) = w(t,l) = . 

Aplica-se o método de separação de variáveis na pesquisa de soluções de (*) : procu ra- 
se soluções da forma 

T(t)X(x) , 
para as quais 

f t T(t)X(x) = t(t)X(x) e 

£,T(t)X(x) = T{t)X^{x). 
Substituindo na equação, e assumindo que T(t)X(x) ^ 0, fica 

(2)/^ _ f{t) _^ 2 X^{X) 



T(t)X(x) = a 2 T{t)X (I >{x) ^^ -±4- = a 



T(t) X(x) 

k k 

onde cada um dos membros tem que ser uma constante real k (porque o membro 

esquerdo não depende de x, o membro direito não depende de t e eles têm que ser 

iguais). 

A solução geral de 

f(t) = jfcT(t) 
é 

T(t) = ce kt com ceK. 

A solução geral de 

X (2 \x) - -\x(x) = 

é, consoante os casos, 

X(x) = c\e^ x + C2e~^ x se k > 

X(x) = ci + c 2 x se k = 

,X(x) = cicos(^^x) + c 2 sin(^^a;) se k < . 
Para que a condição homogénea na fronteira, 

T(t)X{0) = T(t)X(l) = , 
seja satisfeita por uma solução não identicamente nula, tem que ser 

X(0) = X{1) = 

(já que, se X(0) + ou X(l) ^ 0, teria que ser T(t) = 0, \/t). Quando k > 0, a 
única solução X{x) que satisfaz esta condição na fronteira é a solução identicamente 
nula, ou seja, nesses casos, 

X(0) = X(l) = =^> d = c 2 = . 

Quando k < 0, a condição X(0) = impõe c\ = e depois a condição X(l) = 
impõe a uma solução não identicamente nula que 

sin(- ) = . 

a 
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11 



Esta condição pode ser satisfeita por uma solução não identicamente nula desde que 



-k 



= mr 



Q 



n= 1,2,3,... 



ou seja, 

k = —n 2 7r 2 a 2 , n = 1, 2, 3, . . . 

Obtêm-se assim todas as seguintes soluções para o problema homogéneo (*) : 

w n (t, x) = e~ n *■ a t sm(mrx) , n = 1, 2, 3, . . . 



T(í) 



X(x) 



definidas para t > O e para a; e [0,1]. Por linearidade, qualquer combinação linear 
finita desta soluções é ainda solução e, mais geralmente, qualquer série uniformemente 
convergente da forma 

oo 

y^b n w n (t,x) 

re=l 

é a/nda solução. 

Obtêm-se assim as seguintes funções que verificam todas as condições do problema 
posto excepto a condição inicial: 

oo 

20 + 40x + ^2 b n w n (t,x) . 

71=1 

Finalmente, impondo a condição inicial, u(0,x) = 75, terá que ser 

oo 

20 + 40a; + ^ 6,^(0, x) = 75 , 

71=1 

ou seja, 

oo 

y b n sm(nirx) = 55 — 40x . 

71=1 

Para encontrar as constantes b n adequadas, desenvolve-se a função 55 — 40a; em 
série de senos. Isso é equivalente a estender esta função ao intervalo [—1,1] como 
função ímpar, e desenvolver a extensão em série de Fourier. Os coeficientes a n de 
uma função ímpar anulam-se e os coeficientes b n , para n > 1, são dados por 

-í 
b n = I f(x) sin(n7ra;) dx 



= 2 / f(x) sin(n7ra;) dx 

Jo 



= 2 / (55 — 40a;) sin(n-7rx) dx 

'o 



= 2 



cos(n7ra;) , cos(n7ra;) , sin(n7ra;) 
-55 ^ + 40a; ^-40- 



nir 



IÍTT 



n 2 7r 2 



= 2 



-55^ '— + 55— + 40^ '- 



mr 



nir 



nn 



110-30(-l)" 

71 7T 
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Conclui-se que a solução formal do problema posto é 



20 + A0x + J2 110 "t i" 1)n e- 2 - Wí sin(nvrx) . 



n=\ 



D 



Comentário: O problema anterior diz respeito a evolução da distribuição de temperatura 
numa barra situada entre os pontos x = e x = L. No instante inicial todos os pontos da 
barra estão ã mesma temperatura (75) e as temperaturas nas extremidades são mantidas 
constantes igual a 20 em x = e a 60 na extremidade x = 1. <> 



(10) 



Determine a solução do seguinte problema para a equação do calor: 



du a d 2 u 

dt ~ ^dx 2 



se < x < 1 
75 se 1 < x < 2 

para í > OeO <z<2 (satisfazendo a equação diferencial para < x < 2). 



u(x, 0) 



d "(0,t) 



Seja c um parâmetro real positivo. Determine a solução do seguinte problema de 
valor na fronteira para a equação das ondas: 



(11) 



d 2 u _ 2d 2 u 
dt 2 dx 2 



u(x, 0) 



XCOS 



%(x,Q) = u(Q,t) = u(l,t) = Q 
para í>0e0<z< 1 (satisfazendo a equação diferencial para < x < 1). 



(12) 



Use o método de separação de variáveis para determinar uma solução do seguinte 
problema para a equação de Laplace: 



a 2 u 

dx 2 



& 

<->!) 



2=0 



^.(x,0) = f y (x,l) = p x (l,y) = 



dy 

du 

dx 



^(0,y) = sin(27ry) 



para < x < 1, < y < 1. 



Método de separação de variáveis 



(13) 



(a) Determine as soluções para í>0e para x e [0,7r] de 



efe 
c>í 



9 2 M 

9s 2 



u(í,0) = u(í,7r) = 

(satisfazendo a equação diferencial para x g]0, 7r[). 
(b) Determine a solução que satisfaz a condição inicial 

14 (0, X) = ("7T — X)X . 



Resolução: 

(a) Pe/o método de separação de variáveis, procura-se soluções da EDP da forma 

u(t,x) = T{t)X{x) , 
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para as quais 

gT(t)X(x) = f(t)X(x) e 

&T{t)X{x) = T(t)X^(x). 
Substituindo na equação, e assumindo que T(t)X(x) ^ O, fica 

t(t)X(x) = T(t)X^(x) - T{t)X{x) <=> ^|| = ^^ _ i 

k k 

onde cada um dos membros tem que ser uma constante real k (porque o membro 

esquerdo não depende de x, o membro direito não depende de t e eles têm que ser 

iguais). 

A solução geral de 

f(t) = kT(t) 
é 

T(t) = ce kt com ceK. 

A solução geral de 

X {2) {x)-{k + l)X{x) = 

é, consoante os casos, 

X(x) = c 1 e^+ Ix + c 2 e-^+ Ix se k + 1 > O 

X{x) = ci + C2X se k + 1 = O 

X(x) = ci cos(\/—k — lx) + C2 sin(-\/— k — lx) se k + 1 < O . 

Para que a condição na fronteira, 

T{t)X{0) = T(í)X(tt) = O , 

seja satisfeita por uma solução não identicamente nula, tem que ser 

X(0) = X(tt) = O 

(já que, se X(0) ^ O ou X(ir) ^ O, teria que ser T(t) = O, Vt). Quando k + 1 > O, a 
única solução X{x) que satisfaz esta condição na fronteira é a solução identicamente 
nula, ou seja, nesses casos, 

X(0) = X(tt) = => d = c 2 = O . 

Quando k + 1 < O, a condição X(0) = O impõe c\ = O e depois a condição X(ir) = O 
impõe a uma solução não identicamente nula que 

sin(v / — k — Itt) = O . 

Esta condição pode ser satisfeita por uma solução não identicamente nula desde que 

V—k — Itt = nir , n = 1, 2, 3, . . . 



ou seja, 



k = — 1 — n , ji= 1,2,3,... 



Obtêm-se assim todas as seguintes soluções para o problema de valor na fronteira 
dado: 

u n (t, x) = e" (1+ ^ sin(nx) , n = 1, 2, 3, . . . 
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definidas para í>Oe para x e [0,7r]. Por linearidade, qualquer combinação linear 
finita destas soluções é ainda solução e, mais geralmente, qualquer série 

oo 

y^b n u n (t,x) 

n=\ 

é uma solução formal. 
(b) Para que a solução geral da forma acima satisfaça a condição inicial 



tem que ser 



ou seja, 



u(0,x) = (n — x)x , 



^ b n u n (0, x) = (-k - x)x , 

71=1 



y b n sin(nx) = (ir — x)x . 

Para encontrar as constantes b n adequadas, desenvolve-se a função (ir — x)x em série 
de senos. Isso é equivalente a estender esta função ao intervalo [—ir, ir] como função 
ímpar, e desenvolver a extensão em série de Fourier. Os coeficientes da série são 

b n = | / (ir — x)xsm(nx) dx 
n J o 

= 2 / x sin(nx) dx / x sin(na;) dx 

Jo ttJo 

= 2 [-^cos(nx) + ^sin(ra)]^ 

~ cos(nx) + ^f sin(ra) + \ cos(nx) 



i ^ sen e impar 



ÍF[1 -(-!)"] 

8_ 

3 71 

10 se n é par 

(onde as primitivas foram obtidas por primitivação por partes, uma vez no primeiro 
integral e duas no segundo). Logo, a solução pretendida é 



Y 4-e" (1+íl2)í sm(nx) 



n^ir 

71 = 1 

n ímpar 



8 



= £ ,9 Znã e - (l+(2 " +1)2) ^in((2n + l)x). 
^ (2n + l) d 7r 

n=(J v ' 



D 
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(a) Determine as soluções para í>0e para x e [0, 1] de 

du 9 2 - 



at dx 2 "T" u 



(14) 



u(t, 0) = , u(t, 1) = sin 1 

(satisfazendo a equação diferencial para x e]0, 1[). 
(b) Determine a solução que satisfaz a condição inicial 

ít(0, x) = 3 sm(27nr) — 7 sin(47nr) + sin a; 



Resolução: 

(a) /A solução geral deste problema com condições não homogéneas pode ser obtida 
somando a uma solução particular do problema não homogéneo a solução geral do 
problema homogéneo associado. 

Vai-se procurar uma solução particular entre as soluções estacionárias, i.e., que não 
dependem do tempo t. Em resumo, procura-se escrever a solução geral na forma 

u(t, x) = v(x) + w(t, x) , 

onde v(x) é uma solução particular estacionária do problema dado e w(t,x) é a 
solução geral do problema homogéneo 

{dw _ d 2 w | 
dt - dx 2 "T" W 
w(t,0) = w(t,l) = . 

Solução particular estacionária: 
Substituindo v(x) na equação, obtém-se 

= t) {2) + v 
cuja solução geral é 

v(x) = c\ cos x + C2 sin x . 

Para que a condição em x = 0, v(0) = 0, seja satisfeita, tem que ser c\ = 0. Depois 
para que a condição em x = 1, v(l) = sinl, seja satisfeita, tem que ser C2 = 1. 
Conclui-se que 

v(x) = sina; 

é uma solução particular (estacionária) do problema não homogéneo. 

Solução geral homogénea: 

Pelo método de separação de variáveis, procura-se soluções do problema homogéneo 

(*) da forma 

w{t,x) = T(t)X(x) , 
para as quais 

§- t T(t)X(x) = t(t)X(x) e 

&T{t)X(x) = T{t)X^{x). 
Substituindo na equação, e assumindo que T(t)X(x) ^ 0, fica 



f{t)X{x) = T{t)X^\x) + T{t)X{x) <=> ^| = :^to 



+ 1 



k k 

onde cada um dos membros tem que ser uma constante real k (porque o membro 
esquerdo não depende de x, o membro direito não depende de t e eles têm que ser 
iguais). 
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A solução geral de 

f(t) = kT(t) 
é 

T(t) = ce kt com ceK. 

A solução geral de 

X { - 2 \x)-{k-l)X{x) = 

é, consoante os casos, 

X{x) = c\e^ =lx + c 2 e~^~ [x se k - 1 > 

X{x) = ci + C2X se k — 1 = 

X(x) = ci cos(\/l — kx) + C2 sin(\/l — kx) se k — 1 < . 
Para que a condição na fronteira, 

T(t)X(0) = T(t)X(l) = , 
seja satisfeita por uma solução não identicamente nula, tem que ser 

X(0) = X(l) = o 

(já que, se X(0) + O ou X(l) ^ O, teria que ser T(t) = O, Vi j. Quando k - 1 > O, a 
única solução X{x) que satisfaz esta condição na fronteira é a solução identicamente 
nula, ou seja, nesses casos, 

X(0) = X(l) = =^> ci = c 2 = O . 

Quando k — l<0,a condição X(0) = O impõe ci = O e depois a condição X(l) = O 
impõe a uma solução não identicamente nula que 



sinfVl - k) = O . 
Esta condição pode ser satisfeita por uma solução não identicamente nula desde que 

vi — k = rnr , n = 1, 2, 3, . . . 

ou seja, 

k = 1 — n 2 7r 2 , n = 1, 2, 3, . . . 

Obtêm-se assim todas as seguintes soluções para o problema homogéneo (•) : 

w n (t, x) = e (1 ~" 27r2) * sm(nirx) , n = 1, 2, 3, . . . 

definidas para t > O e para x e [0,1]. Por linearidade, qualquer combinação linear 
finita desta soluções é ainda solução e, mais geralmente, qualquer série 

oc 



y^b n w n (t,. 



n=\ 



é uma solução formal. 

OD 

Solução geral: sin x + \^ b n w n (t, x) . 

ra=l 

(b) Para que a solução geral da forma acima satisfaça a condição inicial 

u(0,x) = 3sin(27nc) — 7sin(47ra;) + sinx , 
tem que ser 

oo 

2_] b n w n (0,x) = 3sin(2-7ra;) — 7sin(47ra;) , 

n=l 
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ou seja, 



y^b n sm(nirx) = 3sin(27ra;) — 7sin(4-7ra;) . 



re=l 

As constantes b n adequadas são 

Ò2 = 3 , Ò4 = — 7 , toc/os os outros b n 's são zero . 

(A função 3sin(27rx) — 7sin(47ra;) já é dada em série de Fourier.) 
A solução pretendida é 

sina; + 3e( 1 - 47r2 ) í sin(2ra) - le^- 16 ^ sm(Airx) . 



D 



(a) Calcule a série de Fourier da função f(x) = x(x 2 — 1) para x e [—1,1]- 

(b) Recorrendo ao método de separação de variáveis, determine as soluções para 
í>0e para x e [0, 1] de 



(15) 



Õ 2 i 



1 di. 



dt 2 ^ * dt 



u = 



(Pi 

O.!- 



u(t,0) = u(t,l) = . 
(c) Determine a solução que satisfaz a condição inicial 



u(0, x) = sin(7nr) 



du 



(0, x) = x(x 2 — 1) 



(16) 



(a) Calcule a série de Fourier da função / 

m = 



definida por 



sm 



7TX 



(b) Recorrendo ao método de separação de variáveis, determine as soluções para 
í>0e para x e [0, 1] de 



du d 2 u 

dt dx 2 

g(í,0) = g(í,l) = 0. 



(c) Determine a solução que satisfaz a condição inicial u(0,x) = sin 



TTX ^ 

2 1 
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Transformada de Laplace 



(17) 



Útil 


zando a 


transformada 


de 


Laplace, 

íy (2) - 
\y(0) 


resolva o seguinte 

- 5y + 4y = e 2t 
= 1, 3/(0) = -1 


problema 


de 


valor 


inicial: 



Resolução: Seja Y(s) a transformada de Laplace da solução y(t). Aplicando a transfor- 
mada de Laplace, a equação fica 

C{yW} - 5£{y} + 4Y(s) = £{e 2t } 
<=> (s 2 Y(s) - sy(0) - y(0)) - 5 (sY(s) - y(0)) + AY{s) = ^ 
^^ s 2 Y(s)-5sY(s) + 4Y(s)-s + 6= ^ 
^^ (s 2 -5s + 4)Y(s) = ^ + s-6 



Y(s) = 



+ 



s-6 



(s-2)(s 2 -5s+4) ^ s 2 -5s+4 " 

Para reconhecer esta função como a transformada de Laplace de uma combinação de 
funções elementares, decompõe-se em fracções simples. 
Cálculos auxiliares : Quanto a primeira parcela de Y(s), tem-se 

s 2 - 5s + 4 = <S=> s = 1 ou s = 4 



( s -l)(s-2)(s-4) 



A | _B_ | _C_ 
s-1 ~r s-2 ~r s-4 



1 = 4(s - 2)(s - 4) + £(s - l)(s - 4) + C(s - l)(s - 2) . 

Quando s = 1, obtém-se que A = |; quando s = 2, obtém-se que B = —\; quando s = 4, 
obtém-se que C = g. /.ogo, 

1 i 

-L Q 



+ 



(s-2)(s 2 -5s + 4) s-1 s-2 s-4' 
Quanto a segunda parcela de Y(s), tem-se 



s-6 



s-4-2 



1 



s 2 -5s+4 (s-l)(s— 4) s-1 (s-l)(s-4) 

1 _ _D_ , _E_ 

(s-l)(s-4) s-1 "T" s-4 



1 = D(s - 4) + £(s - 1) • 
Quando s = 1, obtém-se que D = — |; quando s = 4, obtém-se que E = ^. Assim, 

s-6 1 



+ 



s 2 - 5s + 4 s-1 s-1 s-4 
Continuação da resolução: De acordo com os cálculos auxiliares, 

i i 



y(s) = 



111 2 2 

3 _ 2 _i_ 6 _i_ 1 i 3 _ 3 

S-1 S-2 ~r S-4 "I" S-1 ~r S-1 S-4 



,11 

^ _ 2 2 

s-1 s-2 s-4 



= 2£{e í } - \C{e 2t ) - ±£{e 4í } 
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pelo que 



y [t) = 2e l - -e 2t - -e u . 
ww 2 2 



D 





Sabendo que y = y(t) satisfaz a equação diferencial 




y^-2y + y = f(t) 


(18) 


e que a transformada de Laplace de y(t) é dada por 

2 1 
Y( ") 1 




K " } (s 2 -i)(s-i) (s-iy ' 




determine / e as condições iniciais de y(t), ou seja, os valores de y(0) e y(0). 



Resolução: 

Solução 1: Aplicando a transformada de Laplace, a equação fica 

C{yW}-2C{v} + Y{ 8 ) = C{f{t)} 



( s 2y (s) _ sy{0) _ y (0) ) _ 2 {sY{s) _ y(0)) + Y (s) = £{/(*)} 
( S 2 - 2s + 1) y(s) - sy(0) - jr(O) = £{/(*)} 



(« " l) 2 ( ( ^-i)(, s -i) + çtijt) = £{/(*)} + 2/(0) + sy(0) 
3 | T + l = £{/(í)} + y(0) + S2/ (0) 

^- = £{/(í)} donde se conclui que 2e _í = /(í) 

1 = 2/(0) 
=y(0) 

porque a transformada de Laplace de uma função não pode incluir termos da forma a + bs. 
Portanto, 

' f{t) = 2e- t 

y(0) = 1 e 

.2/(0) = 

Solução 2: Sabendo que 

2 1 

1» = ,„, TT7- ^ + 



(s 2 -l)(s-l) (s-1) 2 ' 

pode-se calcular y(t) desenvolvendo a expressão da direita em fracções simples. Cálculos 
auxiliares: 

(s 2 -l)(s-l) = (s-l)\s + l) 

2 ABC 

+ 



(S 2 -1)(S-1) S-l (S-1) 2 8+1 

2 = A(s- l)(s + 1) + B(a + 1) + C(s - l) 2 

Quando s = 1, obtém-se que B = 1; quando s = — 1, obtém-se que C = \. 
Derivando a expressão acima, 

= A(s - 1) + A(s + 1) + B + 2C(s - 1) , 
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e avaliando, por exemplo, em s = 1, obtém-se que A = — \. Logo, 



Y(s) = 



._2 I 2 I £ 

s-1 ^ s+1 ^ (s-1) 2 



-\C{e t } + \C{e- t }-2lC{e t } 
-|£{e í } + |£{e- í } + 2£{íe í } 



donde se conclui que 
Como 

os valores em são: 

2/(0) = 



1 , 1 



y(t) = --e t +-e- t + 2te t . 



m = y 



1 1 
-2 + 2 = ° 



_í + 2tê 



, x 3 1 



/4 função f(t) pode ser calculada substituindo y(t) na equação: 

„< 2 >(t) = -** + !. 

pelo que 

yW -2y + y 



-é + -e _í + 2te t , 
2 2 



;£ e * + i c -* + 2tê) 



2 (|e* - |e _t + 2íe*) + 



■^6*+ ie" í + 2íe í 



= 2e" 



D 



(19) 



Utilizar 


ido a 


transformada 


de Laplace, 

íy (2) - 
\y(o) 


resolva o seguinte 

-4y + 4y = f(t) 
= 0, 3/(0) = 2 


prob 


ema de valor 


inicial: 


onde 




















/(*) = 


íte 2í , 

|íe 2í +(í- 


l) e 2(*-D j 


se 
se 


< t < 1 
t > 1 . 





(20) 



Utilizar 


ido a 


transformada 


de Laplace 


, resolva o segi 


linte problema 


de 


valor 


inicial: 








íy (2) + 2y + 


y = 2(t- 


3)H 3 (Í) 














\y(o) = 


2. 


2/(0) = 1 


? 








onde 






H 3 (t) = < 


ío, 

l 1 ' 


se 
se 


< í < 3 
t > 3 . 









